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Abstract 



Recent redshift-distance measurements of Type la supernovae (SNe la) at 

cosmological distances suggest that two-third of the energy density of the 

universe is dominated by dark energy component with an effective negative 

pressure. This dark energy component is described by the equation of state 

^ ' p x = wp x (w > —1). We use gravitational lensing statistics to constrain 
o I 

the equation of state of this dark energy. We use n(A0), image separation 
distribution function of lensed quasars, as a tool to probe w. We find that for 
the observed range of Q m ~ 0.2 — 0.4, w should lie between —0.8 < w < —0.4 
•rH ! in order to have five lensed quasars in a sample of 867 optical quasars. This 
?-j ; limit is highly sensitive to lens and Schechter parameters and evolution of 
galaxies. 
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1 Introduction 



Recent redshift-distance measurements of SNe la at cosmological distances 
suggest that 2/3 of the energy density of the universe could be in the form 
of a dark energy component with an effective negative pressure [|I], @, |3|, @|. 
This causes the universe to accelerate. There is other observational evidence 
which also supports the existence of an unknown component of energy density 
with pressure p = wp. This includes the recent measurements of the angular 
power spectrum which peak at / ~ 200 j5| and the dynamical estimates of the 
matter density Q m = 0.35 ± 0.07 (jOj. Many candidates have been proposed 
for this dark energy component which is characterised by an equation of state 
— I < w = p/ p < 0. 

The first candidate is the cosmological constant characterised by to = 
— 1. In this case the vacuum energy density or the A — term (87rG/9 vac ) is 
independent of time ( for a recent review see [|7j). There are several other 
possibilities for this dark energy: 

1. A time varying A — term| 



2. Rolling scalar fields (quintessence) 

3. Frustrated network of topological defects in which w = —N/3 where N 
is dimension of the defect [[TD| . 



4. X-matter ITD 



There are some independent methods for constraining w. For example 
Perlmutter et al. (1999) obtained w < —0.6 (95% CL) using large-scale 
structure and SNe la in a flat universe 0. Waga & Miceli (1999) used lensing 
statistics and supernovae data to find that w < —0.7 (68% CL) in a flat 
[TZfl . Recently, Lima & Alcaniz (2000) used age measurements of old 



universe 



high redshift galaxies (OHRG) to limit w [|T3j. In flat universe for Q TO = 0.3, 
the ages of OHRG give w < —0.27. By combining the "cosmic concordance" 
method with maximum likelihood estimator, Wang et al. (2000) found that 
the best-fit model lies in the range Q m = 0.33±0.05 with an effective equation 
of state w ~ -0.6 ± 0.07 jTf . 

In this work we model this unknown form of dark energy as X-matter 
which is characterized by its equation of state, p x = w x p x and w x > — 1. 
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We use gravitational lensing as a tool to constrain the w for X-matter. In 
section 2 we describe some formulae (age-redshift relationship and angular 
diameter distance) used in lensing statistics. In section 3 we explain the use 
of n(A6>), the image separation distribution function of lensed quasars, as a 
tool to constrain the cosmic equation of state of dark energy. We summarize 
our results in sec. 4. 



2 Cosmology with dark energy 



We consider spatially flat, homogeneous and isotropic cosmologies. The Ein- 
stein equations in the presence of nonrelativistic matter and dark energy are 
given by: 



H\z) = C-y = Hi 



04§) 3 +^(§) 3(1 - 



a 



^) 3 + (i + 3w) n x 



Rq 

where the dot represents derivative with respect to time. Further 



(i) 



(2) 



where Hq is the Hubble constant at the present epoch, while p m o and p X Q are 
the nonrelativistic matter density and the dark energy density respectively 
at the present epoch. For more details see ref . [|T5| . 



The age of the universe at the redshift z is given by 

dz' 



Hot(z) = H jf°°- 



+ z')H(z') 



oo 



dz' 



(i + z yn m {i + z'f + 0,(1 + z'f^) 

The angular diameter distance between redshifts Z\ and z% reads as 

, / \ f Z2 dz 

dA{Zi,z 2 ) = 



(3) 



i + z2)J* Jn m a + z)s + n x (i + zW+ w ) 



(4) 
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3 n(AO) as a probe 



n(A6) is the image separation distribution function for lensed quasars. To 
understand n(A6) we first have to calculate the optical depth. The lensing 
probability or the optical depth dr of a beam encountering a lensing galaxy 
at redshift zl in traversing dzi is given by the ratio of the differential light 
travel distance cdt to its mean free path between successive encounters with 
galaxies (l/nL(z)a), til{z) is the number density of galaxies and o is the 
effective cross-section for strong lensing events. Therefore, 

ait 

dr = nAz)o- — dz^. (5) 
dz L 

We further assume it to be conserved: til{z) = n (l + ^l) 3 . The present-day 
galaxy luminosity function can be described by the Schechter function |T6j 

$(L, z = 0)dL = ^ (£) Q exp (-£) £. (6) 
The present day comoving number density of galaxies can be calculated 

as 

n = f <S>(L)dL. (7) 

The Singular Isothermal Sphere (SIS) model with one dimensional velocity 
dispersion v is a good approximation to account for the lensing properties of 
a real galaxy. The deflection angle for all impact parameters is given by 
a = 47T-u 2 /c 2 . The lens produces two images if the angular position of the 
source is less than the critical angle (3 cr , which is the deflection of a beam 
passing at any radius through an SIS: 

(3 cr = aD LS /D os , (8) 

We use the notation D l = d A (0, z L ), D LS = d A (z L , z s ), D s = <U(0, z s ). 
Then the critical impact parameter is defined by a cr = DolPct and the cross- 
section is given by 
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The differential probability dr of a lensing event can be written as 



dr 



dz 



L 



16?r 3 



DnrD 



cdt 



" 4 U^r) ^ dT L dZL - (10) 

The total optical depth can be obtained by integrating dr from to zs which 
is equal to [[TY| 



7 — 1 

rfe)= 3Q (D os (l + z s )f(R )- 3 



with 



f* = — -3- 0* vjr + - + 1 



We neglect the contribution of spirals as lenses as their velocity dispersion 
is small as compared to ellipticals. The relationship between the luminosity 
and velocity is given by the Faber-Jackson relationship j- = 

Table 1 lists lens and Schechter parameters as given by Loveday et al. 
[|TBf (hereafter LPEM parameters) and Kochanek |T7j (hereafter K96 param- 
eters) . For most of our calculations we use LPEM parameters unless specified 
otherwise. 

The differential optical depth of lensing in traversing dzi with angular 
separation between (f) and (j) + d(j) is given by 



d 2 r 

dzidcj) 



dz T , = 



xF 



Dos 



ID LS 



D 



exp 

( DqlDls 

\ RqDqs 



OS 



D LS 

2 



1 

Rq 



cdt 
dz L 

dz L 



'12: 



The normalized image angular separation distribution for a source at zs 
is obtained by integrating the above expression over z^. 



dV 



1 



d 2 T 



-dz 



d(f) r(z s ) J o dz L d(f~ L ' ^ 
We include two correction factors in the probability of lensing: (1) the 
magnification bias and (2) the selection function due to finite resolution and 
dynamic range. 
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The magnification bias B is an enhancement of the probability that a 
quasar is lensed. The bias for a quasar at redshift z with apparent magnitude 
m is written as 

B(m, A h A 2 , z) = // 2 ^(m + 2.5 log A, z) P (A)dA (14) 

am J M dm 

where p(A) is the probability distribution for a greater amplification A which 
is 8 /A 3 for SIS model. We use A\ = 2 and A^ = 10 4 . We use the quasar 
luminosity function as given by Kochanek |I7(| , 



e/m ^ ' 



where 



m (2) = 



m + {z + 1) if z < 1 

m if 1 < 2 < 3 (16) 

m - 0.7(2 - 3) if z > 3. 

We use x = 1.07, y = 0.27 and mo = 18.92. We considered a total of 862 
'z > 1) highly luminous optical quasars plus five lenses (jT2J . 



Selection effects are caused by limitations on dynamic range, limitations on 
resolution and presence of confusing sources such as stars. Therefore we must 
include a selection function to correct the probabilities. In the SIS model the 
selection function is modeled by the maximum magnitude difference Am(6) 
that can be detected for two images separated by A6. This is equivalent to a 
limit on the flux ratio (/ > 1) between two images / = lo°- 4Am w. The total 
magnification of images becomes Af = Aq(J + l)/(f — 1). So the survey can 
only detect lenses with magnifications larger than Af. This sets a lower limit 
on the magnification. Therefore, A\ in the bias function gets replaced by 
Af(0). To get selection function corrected probabilities we divide our sample 
into two parts: the ground based surveys and the HST Snapshot survey. 
We use the selection function as suggested by Kochanek |j20(| . The corrected 



image separation distribution function for a single source at redshift zs is 
given as [0, ^TJ 
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P(A6) = B(m,A f (6),A 2 ,z)^J Q ZS 



'Dos, 




(Dos\*~ 


[D LS . 


exp 


[ KdJ) \ 



x F* 



cdt (l + z L f 
d *L r ( a + ^ + 1 



DqlDls* 
. RqDqs - 



1 

i?0 



(17) 



Similarly the corrected lensing probability for a given source at redshift z is 
given as [0, gTj 



P 



dP 



B(m,A f (0),A 2 ,z) 



Here and A# are linked through = 



A0 



8ir(v*/c) 2 - 

The expected number of lensed quasars is 



(18) 



T,Pi, where Pi is the 



lensing probability of the ith quasar and the sum is over the entire quasar 
sample. Similarly, the image-separation distribution function for the adopted 
quasar sample is n(A9) = J2Pi{A6). The summation is over all quasars in a 
given sample. 



4 Results and Discussions 

Gravitational lensing statistics is a sensitive cosmological probe for determi- 
nation of the nature of dark energy. This is because statistics of multiply 
imaged lensed quasars can probe the universe to a redshift z ~ 1 or even 
higher. This is the time when dark energy starts playing a dominant role in 
the dynamics of universe. 

That lensing statistics can be used as a tool to constrain various dark 
energy candidates has been known for some time. Kochanek (1996) gave a 
2cr upper bound on 11a < 0.66 from multiple images of lensed quasars [jT7| . 



Waga & Miceli (1999) used the combined analysis of gravitational lensing 
and Type la supernovae to constrain the time dependent cosmological term 
[ |T2l . Their combined analysis shows that w < —0.7. Cooray and Huterer 
(1999) also used lensing statistics to constrain various quintessence models 
[|22fl . However, there are several uncertainities involved in using gravitational 
lensing statistics tool to probe cosmology |T|j . 



In this article we neglect the role of spirals as lenses as their velocity 
dispersion is small as compared to E/SO galaxies. The constraints obtained 
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from lensing statistics strongly depend on lens parameters. There is a need for 
updated Schechter parameters, the faint end slope a and normalisation 0* for 
E/SO galaxies. Earlier work on lensing statistics used a ~ — 1, which implies 
the existence of numerous faint E/SO galaxies acting as lenses. Because of 
limited resolution, this faint part of the luminosity function is still uncertain. 
Also the parameters should be determined in a highly correlated manner from 
a galaxy survey and use of parameters derived from various surveys might 
introduce error. We use the updated luminosity function of LPEM. The 
LPEM luminosity function is characterised by the shallow slope a at faint 
end and the smaller normalisation 0* which shifts the distribution to large 
image separations fl21fl . Fig. 1 shows the expected number of lensed quasars 
n(A9 < 4") as a function of Q m in a flat universe with w = — 1 for the LPEM 
and K96 [jTTj lens and Schechter parameters (see Table 1). We observe K96 
parameters predict more lensed quasars for this sample and for no value of 
the parameter Q m the expected number of lensed quasars becomes equal to 
5. Moreover, as pointed out by Chiba et al. []2Tj] K96 parameters have been 
derived from various galaxy surveys and hence lack consistency. 

Recently several galaxy surveys have come up with a much larger sample of 
galaxies. This has improved our knowledge of the galaxy luminosity function. 
But these surveys don't classify the galaxies by their morphological type [|23[ 



In lensing statistics we need Schechter parameters of early type galaxies only. 
We feel at the moment LPEM survey provides the most complete information 
regarding the early type galaxies. 

We use the image separation distribution function function n(A6) to con- 
strain the cosmic equation of state for the dark energy. n(A6) depends upon 
w through the angular diameter distances as shown in Section 2. By varying 
w, the distribution function changes which on comparison with the observa- 
tions gives a constraint on w. 

Fig.l shows the expected number of lensed quasars n(A6 < 4") as a func- 
tion of Q m in a flat universe with w = — 1. Comparing the predicted numbers 
with the observed lenses, we find a value of Q m = 0.45. We further generate 
10 4 data sets (quasar sample) using the bootstrap method. We find 'best fit' 
Q m for each set to obtain error bars on Q m . We finally obtain Q m = 0.5 ±0.2. 
Therefore the gravitational lensing statistics do not favour a large cosmolog- 
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ical constant. 

In Fig. 2 n(A9) is plotted against A9 (image separation) in the flat cos- 
mology for various values of w. The plotted rectangles indicate the image- 
separation distribution of the five lensed quasars in the optical sample con- 
sidered in this calculation. As indicated by recent distance measurements of 
Type la supernovae, we fix Q m = 0.3 in Fig. 2. On comparing the theoretical 
prediction of image distribution function with the observations we see that 
w = — 1 predicts a large number of lenses which is not supported by the ob- 
servations. On the other hand, w = gives too few lensed quasars. We also 
plot the expected image seperation distribution for w = —0.7 and w = —0.5 
(which corresponds to point (-0.5,0.3) on the curve in Fig. 3). Unfortunately, 
we cannot do a statistical study here as the number of lensed quasars in this 
sample is only five. 

If we increase the value of Q m in a flat universe, a smaller value of w is 
required to match with observations. The magnitude of the peak n(A6) is 
sensitive to the value of Q m in a flat universe. A larger value of Q m for a fixed 
value of z gives a smaller angular diameter distance and hence a smaller value 
of the peak. The position of the peak of n(A6) is sensitive to the value of a 
i.e. the faint end slope of the luminosity function. If we take the conventional 
value of a = — 1, the peak will shift to a lower value of A6 or in other words 
it will predict lensed images with smaller angular separations |Z4j. 

In Fig. 3 we show the Q m -w plane. Each point on the curve corresponds 
to a pair of (Q m , w) for which the expected number of lensed quasar is equal 
to five (til = 5). If the matter contribution Q TO increases, a smaller value of w 
is required to produce five lenses. The dotted lines in this figure correspond 
to the observed range Q m ~ 0.2 — 0.4 |j25fl. The corresponding range for w 



in the flat universe is —0.8 < w < —0.4. These results agree very well with 
the limit obtained on w by other independent methods as described in the 
introduction. With this range of w we can rule out atleast two dark energy 
candidates. First, cosmic strings (w ~ —0.33) and second the cosmological 
constant (w = —1.0). It is quite interesting to compare this result with 
the limit obtained on w by using the MAXIMA-1 AND BOOMERANG- 
98 data. This data treats the dark energy as a quintessence giving a limit 
— l<w< —0.5. Our constraint on w is much tighter than that obtained by 
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using the MAXIMA-1 AND BOOMERANG-98 j26j. 

The constraints obtained here depend strongly on lens and Schechter pa- 
rameters, evolution of galaxies and, of course, on the quality of the lensing 
data. The dependence of n(A6) on lens and Schechter parameters (as shown 
in eq.(12) comes mainly through F* (eq.ll), exponential dependence on 7 
and other related factors present in the expression for the differential opti- 
cal depth in traversing distance dz^ and with angular separation d(j). The 
presence of galaxy evolution decreases the optical depth and hence the con- 
straint on w becomes weaker []27j. The main difficulty is that the number of 



observed lensed QSOs in this sample is too small to put strong constraints on 
w. Extended surveys are required to establish n(A9) as a powerful tool. The 
upcoming Sloan Digital Sky Survey which is going down to la magnitude 
limit of ~ 23 will definitely increase our understanding of lensing phenomena 
and cosmological parameters. Moreover, a large number of new gravitational 
lens systems (~ 18) have been discovered by the Cosmic Lens All-Sky Survey 
(CLASS) [291 . CLASS is the largest survey of its kind which has more 



than 10000 radio sources down to a 5 GHz flux density of 30 mJy But the 
major disadvantage of a radio lens survey is that there is little information on 
the redshift-dependent number-magnitude relation. This may lead to serious 
systematic uncertainties in the derived cosmological constraints. 

We need a larger and complete sample of high-z QSOs, a better under- 
standing of the formation and evolution of galaxies over wide range of red- 
shifts and accurate luminosity function parameters before any definitive and 
strong statements can be made regarding the constraints on w by gravita- 
tional lensing. 
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Figure 1: Predicted total number of lenses with A9 < 4" as function of Vt m for flat universe 
with w = — 1 (constant A). The observed five lensed QSOs from optical lens survey is shown 
by horizantal line. 
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Figure 2: The expected distribution of lens image separations with Vt m = 0.3, Q x = 
compared with the observed image-separation distribution in the optical sample (histogrs 
The expected distribution is plotted for different values of w. 
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Figure 3: Contour between Q m and w in order to get five lensed quasars in the optical 
sample. The dotted lines indicate the observed range of il m ~ 0.2 — 0.4. The vertical lines 
give the corresponding value of w. 
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